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ABSTRACT 



00 

^^ ' We have mapped the stellar and gaseous kinematics, as well as the emission-line flux 

f) I distributions and ratios, from the inner a;450pc radius of the Seyfert 2 galaxy Mrk 1157, 

• , using two-dimensional (2D) near-IR J- and Ki-hand spectra obtained with the Gemini NIFS 

"^ ■ instrument at a spatial resolution of «;35 pc and velocity resolution of k 40 km s ' . 

The stellar velocity field shows a rotation pattern, with a discrete S-shaped zero velocity 
Q . curve - a signature of a nuclear bar The presence of a bar is also supported by the resid- 

H ' ual map between the observed rotation field and a model of circular orbits in a Plummer 

c/3 I potential. The stellar velocity dispersion (cTt) map presents a partial ring of low-cr, values 

^^ ' (50 - 60kms"') at 250 pc from the nucleus surrounded by higher cr, values from the galaxy 

bulge. We propose that this ring has origin in kinematically colder regions with recent star 
pg , formation. The velocity dispersion of the bulge (lOOkms"') implies in a black hole mass of 

>■ Mbh ^^3+lJx 10^ Mq. 

Tlj" I Emission-line flux distributions are most extended along PA= 27/153°, reaching at least 

^O i 450 pc from the nucleus and following the orientation observed in previous optical emission- 

•/^ ' line [O III] imaging and radio jet. The molecular H2 gas has an excitation temperature Tgxc ~ 

Cn I 2300 K and its emission is dominated by thermal processes, mainly due to X-ray heating 

[>--. ■ by the active nucleus, with a possible small contribution from shocks produced by the radio 

f^ ' jet. The [Feii] excitation has a larger contribution from shocks produced by the radio jet, as 

evidenced by the line-ratio maps and velocity dispersion map, which show spatial correlation 
with the radio structures. The coronal lines are resolved, extending up to x 150pc and are also 
slightly more extended along PA= 27/153°. 

The gaseous kinematics shows two components, one due to gas located in the galaxy 
^^ ' plane, in similar rotation to that of the stars and another in outflow, which is oriented close 

H I to the plane of the sky, thus extending to high latitudes, as the galaxy plane is inclined by 

^45° relative to the plane of the sky. The gas rotating in the plane dominates the H2 and 
Pa/3 emission, while the gas in outflow is observed predominantly in [Fe 11] emission. The 
[Fe n] emission is originated in gas being pushed by the radio jet, which destroys dust grains 
releasing the Fe. From the outflow velocities and implied geometry, we estimate an outflow 
mass rate of Mout ~ 6 Moyr"' for the ionised gas and a kinetic power for the outflow of 
Ek23x lO'*' ergs-i a; 0.15 x Lboi- 

The distinct flux distributions and kinematics of the H2 and [Fe 11] emitting gas, with the 
former more restricted to the plane of the galaxy, and the later tracing the outflows related 
to radio jets is a common characteristic of the 6 Seyfert galaxies (ES0428-G14, NGC4051, 
NGC 7582, NGC 415 1, Mrk 1066 and now Mrk 1 157) we have studied so far using similar 2D 
observations, and other 2 (Circinus and NGC 21 10) using long-slit observations. We conclude 
that the H2 emission surrounding the nucleus in the galaxy plane is a tracer of the gas feeding 
to the active nucleus while the [Fe 11] emission is a tracer of its feedback. 

Key words: galaxies: individual (Mrk 1157) - galaxies: individual (NGC 591) - galaxies: 
Seyfert - galaxies: ISM - infrared: galaxies - galaxies: kinematics and dynamics 
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1 INTRODUCTION 

Recent resolved imaging and spectroscopic studies of the cen- 
tral regions of nearby active galaxies have been allowing us 
to probe the feeding and feedback mechanisms of their central 
engines. The study of the of the ionized gas of the Narrow- 
Line Region (NLR) of active galaxies allows to investigate how 
the radiation and mass outflows from the nucleus interact with 
the circumnuclear gas, affecting its kinematics and e xc itation 



(e.g. iFischer et al." 2011, 201 0l: Icrensha w et al.' ^OlOb'b. 200S; 



Crenshaw & Kr aemer 2007; Krae mer et al..2009.:.Holt et al 2006 : 
Veilleux. Cecil & Bland-Hawthorn 2005l:ISchmitt & Kinnevli99e ; 



Veilleux. Goodrich & Wilson 1 1997i : IWilson et. al.l Il993h . On the 



other hand, the study of molecular and low ionization gas can 
tell us about the feeding of the Active Galactic Nucleus (AGN) 
of these galaxies by the mapping of in flows and estimating mass 
inflow rates towards the centre (e.g. van de Ven & Fathil 201CI: 
Sanchez et al.' '2009'; 'Storchi -Bergmann et alj |2007| ; iFathi et al.l 
.2006; Mundell & Shone 1993). However, most of the above 
studies are based on opti cal observations, which are af- 
fected by dust obscuration jFerruit. Wilson & Mulchaevl I200CI ; 



ovi 



iMulchaev. Wilson & Tsvetanovll996al lbh. a problem that can be al- 
leviated by using infrared lines to map the NLR emission. 

Since 2006, we have been observing the inner few hun- 
dred parsecs of nearby active galaxies with Integral Field Spec- 
trographs at the Gemini Observatory, with the goal of map- 
ping both inflows and outflows around nearby active galactic 
nuclei and try to constrain the mass flow rates. In the near- 
infrared (hereafter near-IR) our main findings have been that 
the molecular (H2) and ionised gases present distinct flux dis- 
tributions and kinematics. The H2 emission gas is usually re- 
stricted to the plane of the galaxy, while the ionised gas extends 
also to high latitudes and is associated with the radio emission 
I Riffel et al ."2006', '2008', '2009"; 'Riff'el, Storchi-Bergmann & NagaJ 
201(1; JRiffel & Storchi-Bergmann 20 1 1 ; Storchi-Bergmann et al. 
2OO9II2OI O). The H2 kinematics is usually dominated by rotation 



in the plane, including in some cases streaming motions towards 
the nucleus, while the kinematics of the ionised gas, and in particu- 
lar of the [Fen] emitting gas, shows, in addition, a strong outflow- 
ing component associated with radio jets from the Active Galac- 
tic Nucleus (AGN). A previou s study using long-slit near-IR data 
dStorchi-Bergmann et alj 19991) of the Seyfert galaxies Circinus and 
NGC 2110, showing higher velocity dispersion and more disturbed 
kinematics for the [Fen] emission as compared with the H2 emis- 
sion also supports this scenario. These results so far suggest that 
the molecular gas can be considered a tracer of the feeding of the 
AGN and the ionised gas a tracer of its feedback. Nevertheless, 
our 2D observations - which provide a complete coverage of the 
kinematics - comprise so far only half a dozen galaxies, and, in or- 
der to compensate for uncertain projection and filling factors, more 
objects need to be studied in detail, in order to allow also the quan- 
tification of the mass inflow and outflow rates. 

In this work, we present the gaseous flux distribution and kine- 
matics of the inner 450 pc of another active galaxy, the Seyfert 2 
galaxy Mrkll57 (NGC 591), for which we were able to mea- 
sure also the stellar kinematics, allowing its comparison with 
the gas kinematics, which can be thus better constrained. This 
is only the second case in which w e could measure the stellar 
kinem atics, the first being Mrk 1066 JRiffel & Storchi-BergmannI 
I2OIII) . Mrk 1157 was selected for this stu dy because: (i) it 
presents strong near-IR emission lines ( e.g iNagar et alj 1 19991 ; 
iRiffel. Rodriguez-Ardila. & Pastorizal2006n . allowing the mapping 



of the gaseous distribution and kinematics; (ii) it has extended radio 
emission, allowing the investigation of the role of the radio jet; (iii) 
the K-band CO absorption band hea ds have been observed in pre- 
vious near-IR spectra ( Riffel, Rodriguez-Ardila. & Pa storiza 20()g) 
allowing the measurement of the stellar kinematics. 

Mrk 1157 is an early-type barred spiral galaxy (SBO/a), lo- 
cated at a distance d = 61.1 Mpc, for which 1" corresponds 
to 296 pc at the galaxy. Radio-continuum images at 3.6 cm and 
6 cm show a radio double with extension of I'.'l oriented along 
the po sition angle PA=153° JNagar et alj 1 9991 ; luJvestad & WilsorJ 
Il989r) . Optical images show that the [Oin] emission extends 
up to 4" following the orientation of the radio jet, while the 
Hg emission extends up to 20 " alo ng the east-west direction 
( IMulchaev. Wilson & Tsvetanovll996ai) . The near-IR nuclear spec- 
trum of Mrk 1157 presents strong emission lines of [S m]. He i, H i, 
[Fen] and H2 and the presence of high ionisation species as [Si vi], 
[Svm] and [Six] (e.g. Riffel, Rodriguez-Ardila, & Pastoriza 200g; 
IVeilleux. Goodrich & Wilsonlll997h . Optical and near-IR spectra 
show no evidence of broa d components in the profiles of per- 
mitte d emission lines (e.g. iRiffel, Rodriguez-Ardila. & Pastorizal 



l2006l ; l\^lleux. Goodrich &WilsoiJll997h . but spectropolarimet- 
ric observations reveal the presence of this component in Ha and 
Hfi dMoran et alj200(]h . The near-IR nuclear spectrum of Mrk 1 157 
also shows absorption lines, with the CO absorptions in the H- and 
K-bands being the most prominent ones and including the detec- 
tion of the CN absorption at l.lfxra - a signature of the presence of 
intermediate age stars (Riffel. Rodriguez-Ardila. & Pastoriza 20061 ; 
IRiffel et al.ll2007h . 

This paper is organised as follows. In Sec.|2]we describe the 
observations and data reduction procedures. The results are pre- 
sented in Sec.[3]and discussed in Sec.|4] We present our conclusions 
in Sec.|5] 



2 OBSERVATIONS AND DATA REDUCTION 

Mrk 1 157 was observed with Gemini NIFS dMcGregor et alj2003h 
operating with the Gemini North Adaptive Optics system ALTAIR 
in September/October 2009 under the programme GN-2009B-Q- 
27, following the standard Object-Sky-Sky-Object dither sequence, 
with off-source sky positions since the target is extended, and indi- 
vidual exposure times of 550 s. 

Two sets of observations with six on-source individual ex- 
posures were obtained at different spectral ranges: the first in the 
I-band, centred at 1.25 yum and covering the spectral region from 
1.14yum to 1.36yum, and the second in the Ki-band, centred at 
2.3 yum covering the spectral range from 2. lOyum to 2.53 yum. In the 
I-band, the I_G5603 grating and ZI_G0601 filter were used, result- 
ing in a spectral resolution of » 1.8 A, as obtained from the mea- 
surement of the full width at half maximum (FWHM) of arc lamp 
lines. The Kj-band observations were obtained using the K1_G5607 
grating and HK_G0603 filter and resulted in a spectral resolution 
of FWHM« 3.5 A. In velocity space, the resolution of the observa- 
tions is » 45km s"' for the Kpband and 35 km s"' for the I-band. 

The data reduction was accomplished using tasks contained 
in the nifs package which is part of Gemini iraf package, as well 
as generic iraf tasks. The reduction procedure included trimming 
of the images, flat-fielding, sky subtraction, wavelength and s- 
distortion calibrations. We have also removed the telluric bands 
and flux calibrated the frames by interpolating a black body func- 
tion to the spectrum of the telluric standard star. The final IFU 
data cube in each band contains ~ 4500 spectra, each spec- 
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Figure 1. Top-left panel: HST WFPC2 continuum image of Mrk 1157 obtained through the filter F606W jMalkan. Goriian & Tamil 1998h . Top-right panel: 
2.17/im continuum image obtained from the NIFS data cube. The grey lines are flux contours from each image and the color bars show the flux scale in ai'bitrary 
units for the HST image and in logarithmic units for the NIFS image and the green contours are from the 3.6 cm radio-continuum image of Nagar et al. ( 1999J). 
Bottom panels show typical spectra obtained within an 0'.'25x(y.'25 aperture for the nucleus and for a location at (X.'4 north-west from it (position A). The flux 
is shown in 10"'^ergs"'cm"-A"' units. The box in the HST image shows the NIFS field of view. 



trum corresponding to an angular coverage of (X/OSxC/OS, which 
translates into ~15xl5pc^ at the galaxy and covering the inner 
3"x3"(~900x900pc2) of the galaxy. 

The angular resolution obtained from the FWHM of the spa- 
tial profile of the telluric standard star is 0'.'11±0'.'02 for the J-band 
and a'12±a'02 for the Ki-band, corresponding to 32.6+5.9 and 
35.5±5.9 pc at the galaxy, respectively. 



3 RESULTS 

In the top-left panel of Figure[T] we present an optical im- 
age of Mrk 1157 obtained with the Hubble Space Tele- 
scope (HST) Wide Field Planetary Camera 2 (WFPC2) 
through the filter F606W - effec t ive w avelength/width 
5843 A/1578.7 A JMalkan. Goriian & Tainl Il998h . In the top- 
right panel we present an image obtained from the NIFS data 
cube for the continuum emission around 2.17//m. In the bottom 
panels we present two characteristic IFU spectra integrated within 
a 0'.'25x0'.'25 aperture: the nuclear one and a spectrum from 0'.'4 
north-west of the nucleus (Position A), at the location of the 



north-west component o f the radio double dNagar et al.1 1 19991 : 
IUlvestad&Wilsonlll989l) . 



3.1 Stellar Kinematics 

In order to obtain the stellar line-of-sight velocity distributions 
(LOSVD) we fitted the '-CO and ''CO stellar absorption band 
heads around 2.3 yum in the A "-band spectra using the penal ised 
Pixel-Fitting (ppxf) method of Cappellari & EmsellemI ( |2004|) . as 
explained in iRiffel & Storchi-BergmannI ( 1201 ih . As stellar tem- 
plate spectra we used those of the Gemini library of late spectral 
type stars observed with the Gemini Near-Infrared Spectrograph 
(GNIRS) IFU and NIFS jWinge. RifTel & Storchi-Bergamanr] 
12009^ . The resulting maps for the radial velocity (V,), stellar ve- 
locity dispersion (cr,), and higher order Gauss-Hermite moments 
(/i3, and lu,) are shown in Figure |2] At locations close to the bor- 
ders of the field of view, the S/N ratio of the galaxy's spectra was 
not high enough to allow reliable measurements, which have thus 
been masked out. 

In order to estimate the uncertainties on the measurements of 
v., cr,, hi, and h4, we performed 500 iterations Monte Carlo simu- 
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lations of the kinematical extractions with pPXF,using the Gemini 
Ubrary spectra as templates, as described in lCappellari & EmsellemI 
1120040 ■ Our NIFS spectra have S/N ratio in the 2.3 yum continuum 
ranging from 15 at the borders of the NIFS field to 90 around the 
nucleus. We have performed two Monte Carlo simulations, one 
for S/N = 20 to obtain the uncertainties near the borders of the 
field and another for S/N = 60 to represent the uncertainties at re- 
gions near the nucleus. The output from the Monte Carlo simu- 
lations show that for a S/N = 60 the maximum uncertainties are 
^V, ^ 5kms-', Act. « Skms"', A/13. a 0.04 and hh^. « 0.04, 
while for S/N = 20 they are AV. ~ ISkms"', Act. a 20kms-', 
hJi-i, a 0.05 and Mi^. a 0.04. 

The top left panel of Fig.[2] shows the stellar velocity field, 
from which we subtracted the heliocentric systemic velocity of 
Vj = 4473 ± 8kms"', obtained from the modelling of the stellar 
velocity field (see Sec. l4.lt . The velocity field shows a rotation pat- 
tern with blueshifts to the south-east and redshifts to the north-west, 
with a maximum velocity of a;130kms"' and kinematical centre 
coincident with the position of the continuum peak flux, within the 
uncertainties. 

In the top-right panel of Fig. |2] we present the stellar veloc- 
ity dispersion (cr,) map, which presents values ranging from 30 
to 150 kms"'. A partial ring of low cr, values (» 50kms"') is 
observed surrounding the nucleus at 0'.'6 (~ 180 pc) from it. The 
bottom panels show the higher order Gauss-Hermite moments h^, 
(left) and h^, (right), which measure asymmetric and symmetric 
deviations, respectively, from a Gaussian velocity distribution. The 
values hif and /i4, are small, varying from -0.12 to 0.12, and are 
similar to those observed for ot her galaxies JEmsellem et al.ll2004l : 
iGanda et al.ll2006l . iRifTel et alj 2008, 200i), indicating that there 
are only small deviations of the stellar LOSVD from a Gaussian 
velocity distribution. Nevertheless, some systematic deviations are 
observed in both maps. The /13, seems to be anti-correlated with the 
stellar velocity field, presenting red wings (positive values) mostly 
in the blueshifted side of the galaxy, and blue wings (negative val- 
ues) in the redshifted side. The h^, map has the highest positive 
values (meaning the the LOSVD is more "pointy" than a Gaussian) 
where the red wings are observed and the highest negative values 
(meaning that the LOSVD is less "pointy" than a Gaussian) where 
the blue wings are observed. 

3.2 Emission-Line Flux Distributions and Ratios 

We measured the fluxes of 29 emission lines from [P 11], [Feii], He i. 
He n, H I, H2, [S ix], [Si vn] and [Ca vm], which are listed in Table[T] 
for the nucleus and position A. The strongest emission lines are 
identified in the nuclear spectrum of Fig.[T] In the K-band nuclear 
spectrum we have also identified the CO stellar absorption band 

heads around 2.3 pm. 

We used our routine profft JRitfell I2OI0I) to fit the pro- 
files of [Pn]^l. 1886pm, [Fen] ^ 1. 2570pm, Pa^, H2^2.1218pm, 
[Caviii] /12.3220 pm and [Sivii] /I2.4833pm emission lines with 
Gauss-Hermite series. We have integrated the flux under the 
emission-line profiles after subtracting the underlying continuum 
(obtained from two spectral windows, one at each side of the pro- 
file) in order to map of the flux distributions. These particular lines 
have been chosen because they have the highest signal-to-noise 
(S/N) ratios among their species (coronal lines, forbidden and per- 
mitted ionised gas lines and molecular lines). As observed in the 
nuclear spectrum of Fig. [T] the [Caviii] emission line is partially 
superimposed on the '^CO 3-1 /i2. 322 pm band head. In order to 
correct for the effect of the absorption, we measured the flux of 



[Cavm] emission line after the subtraction of the contribution of 
the underlying stellar population, obtained from the fitting of the 
stellar kinematics (see Sec. 13. 11 1. 

In Figure [3] we present the resulting flux distributions, which 
have typical uncertainties smaller than 10%. The central cross 
marks the position of the nucleus, defined as the location of the 
peak of the continuum emission, and the green contours over- 
laid on the [Fen] and Pa/ 3 flux map s are f rom the 3.6 cm radio- 
continuum image from Nagar et alj 1 19991). W e used the 20 cm 
radio-continuum image from lNagaret alj < 1 9991) to align the near- 
IR and radio images, under the assumption that both present their 
peak emission at the same position - the nucleus. The uncertainty 
in the alignment is estimated to be smaller than 0'.'2. 

The flux distributions for the [Pn] and [Fen] emission lines 
are similar. They are well correlated with the 3.6 cm radio- 
continuum emission, being more extended along the position an- 
gle of the radio double PA=153°. The peak flux is observed at 
0'.'4 north-west of the nucleus, approximately coincident with the 
north-west component of the radio double. A secondary peak is 
observed at 0'.'2 south-east of the nucleus, at the same position of 
the other component of the radio double (hereafter we will call the 
radio components NW and SE hot spots). 

The Pa/3 emission is also extended along PA=153°, but shows 
some emission also to the east of the nucleus, being less collimated 
than the [P n] and [Fe 11] flux distributions; its emission peaks at 0'.'2 
south-east of the nucleus, at the position of the SE hot spot. 

The H2 emission is more uniformly distributed, extending up 
to (^.'8 arcsec along all directions, but the highest intensity levels 
are also more elongated along the orientation of the radio double. 

The coronal-line emission (traced by the [Cavm] and [Sivii] 
emission lines) is resolved by our observations (the FWHM of the 
flux distribution for both emission lines is a;0'.'4), peaks at the nu- 
cleus and extends up to 0'.'4 (w 120 pc) from it. The [Si vii] coronal- 
line emission in more extended along PA= 153°, as observed for the 
other emission lines. 



3.2.1 Emission-line Ratios 

In order to map the Narrow-Line Region (NLR) extinction and ex- 
citation mechanisms of the [Fe 11] and H2 emission lines, we con- 
structed the flux ratio maps shown in Figure|4] 

The reddening map was obtained from the Pa/J/Bry line ratio 
as 



E{B-V) = A.lA\og 



5.88 \ 

F Pap I F B,-y I 



(1) 



where Fp^p and Fgry are the fluxes of Pafi and Bry emis- 
sion lines, respectively. We hav e used the reddening law of 
ICardelli. Clavton & MathisI ( 11989 1 and adopted the intrinsic ra- 
tio Fpag/Fsry = 5.88 Corr esponding to case B recombination 
JOsterbrock & Ferlanci2006r) . The resulting E(B-V) map is shown 
in the top-left panel of Fig.|4l where the highest E{B - V) values - 
corresponding to the most heavily reddened regions - are those rep- 
resented by the colours light yellow and white. The highest values, 
of up to 1.8, are observed at the nucleus and to north-east, while 
the lowest values, of about 0.4, are observed at the radio hot spots. 
Typical uncertainties are smaller than 0.2, and the average value is 
E(B -V)~ 0.5. 

In the top-right panel of Fig.|4] we present the 
[Feii]/ll.2570pm/Pa/? ratio map, which can be used 
to investigate the excitation mechanism of [Fen] (e.g. 
iRodriguez-Ardila et al.l2004l ; lRodriguez-Ardila. Riffel & Pastorizal 
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Figure 2. Stellar kinematics obtained from the fitting of the CO band heads. Top left: centroid velocity field; top right: velocity dispersion map; bottom left: 
h3, map and bottom right: h4» map. The colour bars show the range of values for the velocity and velocity dispersion in kms"' (top) and for h3 and h^ 
Gauss-Hermite moments (bottom). 



20051: IS torchi-Beremann et al.' '2009"; 'Riffel et al.! l2006l: 



Riffel. Sto rchi-Bergmann & Nagar 2010; Storchi-Bergm ann et al.l 
200a)- Seyfert galaxies present typical values for this ratio between 
0.6 and 2.0, Starbursts have ratios <0.6 and Low-Ionisation Nuclear 
Emission-line Regions (LINERs) show values higher than 2.0 (e.g. 



Rodriguez - Ardila et al. 2004; Rodriguez-Ardila, Riffel & Past orizal 
20051; [storchi-Bergmann et al. 2009 ; Riffel et al. 200^: 



Riffel. Storchi-Bergmann & Nagaij I2OI0I ; IStorchi-Bergmann et al.l 
20091). The lowest values for Mrk 1157 (~ 0.3) are observed near 
the nucleus (closer than =s0'.'4 or 120 pc from the nucleus), while 
the highest values, of up to 3, are observed at a;0'.'4 north-west of 
the nucleus, at the position of the NW radio hot spot. 

Another line ratio that can be used to investigate the [Fen] 
excitation mechanism is [Feii]/11.2570yum/[Pii]/11.8861yum. Val- 
ues larger than 2 indicate that shocks have passed through the 
gas destroying the dust grains, releasing th e Fe and enhanc- 
ing its abundance and thus emi ssion (e.g. Oliva et al. 2001; 
Storchi-Bergmann et al.l [2009; Riff el. Storchi-B ergmann & Nagar 
20101) . We present this ratio map in the bottom-left panel of Fig.|4] 



The lowest values of » 2 are observed at the nucleus, while the 
highest values of up to 10 are seen at 0'.'4 north-west, at the loca- 
tion of the NW hot spot. 

In the bottom-right panel of Fig.|4] we present 
the H2/i2.1218yum/Br7 ratio map, which is useful to 
investigate the excitation of the H2 molecule (e.g. 



Rodriguez-Ardila et al. 2004; Rodriguez-Ardila, Riffel & Pastoriza 
2005; Storchi-Bergmann et al. 2009; Riffel et al. 2006, 20oJ 
20091 ; iRiffel. Storchi-Bergmann & NagaJ l2010h . This ratio 
presents values ranging from 0.6 to 2 for Seyfert galaxies (ej 



1^ 



: al.l2( 



Rodriguez-Ardila et al.l2004 ; Rodriguez-Ardila. Riffel & Pastoriz; 



7^ 
i9;T 



2005';'Storchi-Bergmann et al.'2009';'Riffel et al.l2006 . l200ai2009 
Riffel, Storchi-Bergmann & Nagar 2010). The lowest values, down 
to » 0.3 are observed to the north-east of the nucleus and the 
highest values of up to 5 are observed mostly along a strip running 
from south-west to north-west of the nucleus. 



3.3 Gas kinematics 



The PROFIT routine ( lRiffell2010f) . which we have used to fit the emis- 
sion line profiles, also outputs the centroid velocity (K), velocity 
dispersion (cr) and higher order Gauss-Hermite moments (h^ and 
/!4), which have been used to map the gas kinematics. 

In Figure [5] we present the velocity fields obtained from the 
centroid wavelength of each emission line, with typical mean un- 
certainty of 8 kms"'. White regions in this figure represent loca- 
tions where the S/N was not high enough to allow the fitting of 
the line profiles. Most velocity fields show mostly blueshifts to 
the east (left in the figures) and redshifts to the west, with the line 
of zero velocity running approximately vertically (north-south) for 
Pa/3 and Ht and from the north-west to the south-east for [Feii]. The 
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Table 1. Measured emission-line fluxes (in units of 10"'^ erg s"' cm"^) for the two positions marked in Fig.[T]within 0'.'25xO'.'25 aperture. 



AyadM'^) 


ID 


Nucleus 


Position A 


1.14713 


[Pn]lZ)3-^Pl 


16.39 ± 1.07 


21.37 ±7.20 


1.16296 


Heii7-5 


19.84 ±4.01 


7.71 ±0.78 


1.18861 


[Pll]^D2-^ Pi 


38.48 ± 3.49 


49.68 ± 8.03 


1.19665 


[Feii]b*Ds,2-a^F3/2 


6.00 ± 0.82 


2.48 ± 0.33 


1.19723 


Hei^Di -3po2 


- 


1.79 ±0.50 


1.29763 


[Fen]a'^Di,2-a''Ds,2 


2.53 ± 0.27 


- 


1.23878 


[Feu]c^Gi/2 -a'^G\ii2 


4.64 ± 0.67 


5.19 ±0.91 


1.24211 


[Fen]fo''Di/2-a*^S5/2 


- 


2.13 ± 1.15 


1.25235 


[Six]3p, -3p2 


34.83 ± 1.57 


4.57 ± 1.17 


1.25702 


[Fen]a''z)7/2-fl*£>9/2 


101.07 ± 7.65 


207.82 ± 19.66 


1.27069 


[Feii]fl4Di/2-a«£>i/2 


- 


18.07 ± 3.23 


1.27912 


[Fen]a''D3/2-a*£l3/2 


7.18 ±1.31 


- 


1.28216 


HiPa/3 


123.27 ± 3.05 


101.69 ±31.30 


1.28607 


[Fen]c2z)5/2-&*^G9/2 


- 


1.84 ±0.29 


1.29462 


[Vcii\a^Df,i2-a^Df,i2 


- 


5.85 ± 0.69 


1.32092 


[¥<iii\ a'^Di 12 -a^Di 12 


32.76 ±3.15 


67.50 ± 8.36 


1.32814 


[Fen]a''D5/2-a''D3/2 


13.99 ±5.80 


8.77 ± 7.4 


2.12183 


H2 1-0 S(l) 


74.57 ± 9.48 


50.46 ± 2.43 


2.15420 


H2l-0S(2) 


3.48 ± 3.60 


2.70 ± 0.96 


2.16612 


HiBry 


46.14 ±6.33 


19.63 ± 1.53 


2.22344 


H2l-0S(0) 


22.32 ± 0.90 


13.02 ± 1.08 


2.24776 


H2 2-1S(1) 


8.41 ± 3.43 


6.27 ± 1.36 


2.32204 


[Caviii]2p0^,-2p0^2 


36.83 ±11.25 


5.12 ±0.99 


2.40847 


H2l-0Q(1) 


56.59 ± 2.69 


33.08 ± 0.78 


2.41367 


H2l-0Q(2) 


14.46 ±6.18 


10.23 ± 1.28 


2.42180 


H2l-0Q(3) 


51.63 ±3.48 


33.77 ± 2.02 


2.43697 


H2l-0Q(4) 


16.21 ±3.88 


9.03 ± 1.59 


2.45485 


H2l-0Q(5) 


32.20 ±3.52 


19.82 ± 4.39 


2.48334 


[Sivn]3pi -3p2 


85.95 ± 5.64 


17.88 ± 1.36 
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Figure 3. Emission-line flux distributions The central cro ss marks the positio n of the nucleus, the grey contours correspond to each flux map, while the green 
contours are from the 3.6 cm radio-continuum image of lNagar et alj 119991) . Colour bars show the range of flux values for each emission line in units of 
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Figure 4. Emission-line ratio maps; the reddening E(B-V) was obtained from Pa/ J/Bry line ratio ( top-left panel). The central cross marks the position of the 
nucleus and the green contours are from the 3.6 cm radio-continuum image of iNagar et alj 1199911 . We do not show the borders of field (the outermost CI.' 5) 
due to the small S/N ratio in the line-ratio maps. 



gas velocity fields differ from the stellar one, for which the line of 
zero velocity runs from the north-east to the south-west. The [Fe ii] 
emitting gas presents higher redshifts than the other emission lines, 
of up to 150 km s"' , observed to south-east of the nucleus, near to 
the edge of the radio structure. Similar redshifts are also observed 
in the [Pii] velocity field, which presents additionally some high 
blueshifts. As observed in the bottom-central and bottom-right pan- 
els, the coronal emitting gas, although just barely extended, seems 
to show similar kinematics to that observed for the rest of the ion- 
ized emitting gas. 



Figure|6]shows the velocity dispersion maps, which have typ- 
ical mean uncertainties of 10 kms"'. The green contours overlaid 
to the [Pen] a map are from the 3.6 cm radio continuum image. 
The [Fe ii] a map shows the highest values of up to a;250 km s"' in 
regions at and around the radio hot spots, while the lowest values 
(» 70 kms"') are observed predominantly to north-east of the nu- 
cleus, away from the radio structure. The [Pii] cr map is similar to 
the [Fe ii] a map, although with k30% lower values. The same is 
observed in the Pa/? cr map but with even lower values (50% lower), 
ranging from »50 km s"' (observed to north-east of the nucleus) to 
w 130 kms"', observed along PA=153° in association with the ra- 
dio structure. The H2 emitting gas presents the lowest a values, 
which are smaller than 60 kms ' at most locations. The coronal 
lines have similar cr values to those observed for the [Pii], with 



median values of » 1 10 km s ' , being smaller than those seen in the 
[Fell]. 

In FigureQwe present the h^ maps, which show values rang- 
ing from -0.4 to 0.4 for most emission lines. The mean uncertain- 
ties in h^ are smaller than 0.03 for all emission lines. The h^ maps 
for the [Pii], [Feii], Pa/3 and H2 emission lines show negative values 
- meaning blue wings - in the profiles to the north-west and posi- 
tive values - meaning red wings - to the south-east, with the high- 
est absolute values observed near the borders of the radio structure. 
The coronal lines of [Caviii] and [Sivii] present mostly negative 
values, indicating the dominance of blue wings in their profiles. 

The h^ maps show values close to zero at most locations and 
are thus not shown. The mean uncertainties are smaller than 0.05 
for all emission lines. The [Feii] and [Pii] show a few high values, 
of up to 0.4, to north-east and to south-west - away from the radio 
structure, meaning that the profiles are more peaked than gaussians 
at these locations. 



3.4 Channel Maps 

Channel maps along the emission line profiles are shown in Fig- 
ures[8l[9]and[T0]for the [Feii], Pa/3, and H2 emitting gas. Each panel 
presents the flux distribution in logarithmic units integrated within 
the velocity bin centred at the velocity shown in the top-left corner 
(relative to the systemic velocity of the galaxy). The central cross 
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Figure 5. Centroid velocity fields for the [Pii]/ll.l886/im (top left), [Fen] /tl. 2570 /im (top central), Pa/3 (top right), H2/12.1218/jm (bottom left), 
[Cavni] /I2.3220yum (bottom central) and [Si vn] /i2.483 3 lira (bottom right) emitting gas. The central cross marks the position of the nucleus and the black 
contours are from the 3.6 cm radio-continuum image of iNagar et alj y999). We do not show the outermost (X'S borders of the NIFS field, since it was not 
possible to fit the emission line profiles at these locations due to the small S/N ratio. The color bars show the range of values for the velocities in km s"' . 



marks the position of the nucleus and the green contours overlaid t o 
some panels are from the 3.6 cm radio image of Naaar et alj (19991). 
We do not show channel maps for [P ii] because they are similar and 
noisier than the [Fen] channel maps. We also do not show channel 
maps for the coronal lines because they do not bring additional in- 
formation to that already contained in the velocity fields shown in 
Figs. [5] and |6] 

In Figure [8] the channel maps along the [Fen] emission-line 
profile show the flux distributions integrated within velocity bins 
of 76 km s~' (corresponding to three spectral pixels). The highest 
blueshifts, which reach -400 kms"', are observed at 0'.'4 north- 
west of the nucleus, approximately coincident with the radio hot 
spot, while the highest redshifts with similar velocities (400 km s"') 
are observed to both sides of the nucleus (at 0'.'4 north-west and 0'.'6 
south-east) along the orientation of the radio structure. A good cor- 
relation between the radio image and the [Feii] flux distribution 
is observed at velocities close to zero, where the [Feii] emission 
shows two peaks, at 0'.'4 north-west and 0'.'2 south-east of the nu- 
cleus, at the location of the two radio hot spots. 

Figure[9]shows the channel maps for the Pa/3 emitting gas for 
the same velocity bin as for [Fen]. The highest blueshifts and red- 
shifts, reaching a;300kms"', are observed to both sides of the nu- 
cleus along the orientation of the radio jet (to north-west and south- 
east of the nucleus). The correlation between the radio and line 
emission observed for the [Fe ii] emitting gas for low velocities is 
not present in the Pa/J channel maps, as the emission at low veloci- 
ties (between » -100 and a 100 km s"') moves from the south-east 
(left and bellow the nucleus) to the north-west (right and above 
the nucleus), approximately following the stellar rotation pattern 



(blueshifts at the south-east to redshifts to the north-west of the nu- 
cleus). At the panel centred at -43 km s"' there is also a structure at 
w r.'6 north-east of the nucleus, still observed at the panel centred 
at 32 km s"' although closer to the nucleus, and at the panel centred 
at 106 kms"' there is another structure at » I'.'S north-west of the 
nucleus. 

The channel maps along the H2 emission-line profile are 
shown in Figure [TO] for a velocity bin of 60 km s"' (corresponding 
to two spectral pixels). These maps have some similarities to those 
of Pa/?, in the sense that at low velocities, the blueshifts are mostly 
observed to the east and south-east and redshifts mostly to the west 
and north-west, indicating rotation. But they reach lower velocities 
than Pa/?, are more uniformly distributed around the nucleus and 
are less correlated to the radio structure. 



4 DISCUSSION 



4.1 Stellar kinematics 



The stellar velocity field of the central region of Mrk 1 157 (top-left 
panel of Fig.[2l(, suggests rotation. In order to obtain the systemic 
velocity, orientation of the line of nodes, the location of the kine- 
matical centre and an estimate for the bulge mass, we have fitted the 
stellar velocities with a model of circular orbits in a plane subject to 
a Plummer gravitation al potential, in which the velocity is given by 
( Barb osa et al. 2006 : Riffel et alj2008l : lRiflFel & Storchi-BergmanrJ 
I2OIII): 
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Figure 6. cr maps for the same emiss i on lin es of Fig[5] Tlie central cross marks the position of the nucleus and the green contours are from the 3.6 cm 
radio-continuum image of lNagar et alj 09991) . We do not show the outermost Cf.'S borders of the NIFS field, since it was not possible to fit the emission line 
profiles at these locations due to small S/N ratio. The color bars show the range of o" values in km s"' . 
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Figure 7 . hj maps the same emission lines of Fig|5] The central cross marks the position of the nucleus and the contours are from the 3.6 cm radio-continuum 
image of lNagar et alJ (1999). We do not show the outermost ff.'S borders of the NIFS field, since it was not possible to fit the emission line profiles at these 
locations due to small S/N ratio. The color bars show the range of /13 values. 
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Figure 9. Same as Fig.[8]for the Pa/3 emission line. 
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Figure 10. Same as Fig.[8]for the H2 emission line for a velocity bin of 60 km s 



Vr = V, + 



R^GM 



sinCOcosC^P - "Po) 



(/?2 + A2)3/2 l_,,^ _ ^^ ^ .in^('.-To) x3/4 



(cos^CP - >Po) + ■ 



(2) 



ccs^d) J 



where R is the projected distance from the nucleus in the plane of 
the sky, Y is the corresponding position angle, M is the mass in- 
side R, G is the Newton's gravitational constant, Vs is the systemic 
velocity, i is the inclination of the disc (/ = for a face on disc), 
To is the position angle of the line of nodes and A is a scale length 
projected in the plane of the sky. 

The equation above contains six free parameters [including 
the location of the kinematical centre (Xq, Yq)], which can be deter- 
mined by fitting the model to the observations. This was done using 
a Levenberg-Marquardt least-squares fitting algorithm, in which 
initial guesses are given for the free parameters. In figure [TT] we 
show the model which gave the best fit (top panel) and the corre- 
sponding residual map, obtained from the subtraction of the mod- 
elled from the observed velocities. The residuals are smaller than 
30 kms"' at most locations, indicating a good reproduction of the 
velocity field, but the residual map shows some systematics. In par- 
ticular, the "blueshifted strip" crossing the middle of the residual 
map and the partial "redshifted strip" to the left of it, resemble the 
residuals we have observed in the stellar velocity field of Mrk 1066, 
whi ch we have attributed to the prese nce of an oval distortion or 
bar ( IRifl^el & Storchi-Bergmannll201 11) . This may also be the case 
of Mrk 1 157, since it seems to have also a nuclear bar, as found by 
iMalkan. Goriian & T am ( 1998) using HST broad band images. The 
stellar velocity field (Fig.HJ also suggests the presence of a small 
S-distortion in the "zero-velocity curve", a signatu re of the pres- 
ence of an asymmetry in the gravitational potential jCombes et al.l 
ll995l ; lEmsellem et alj20o3 : lRiffel & Storchi-Bergmanrj|201 ll) . 



The parameters derived from the fit are: the systemic velocity 
corrected to the heliocentric reference frame Vi =4473 ± 8 km s"', 
To =112° ±5°, M =5.4±0.5 X 10** Mo, i = 46° ± 6°, and 
A =119± 18 pc. The derived kinematical centre is Xq = -1.5±13pc 
and }o=69±13 pc measured relative to the position of the peak of 
the continuum emission. The sy stemic velo city is abou t 25 km s~' 
smaller than the one quoted i n IVeron-Cett v & VerorJ ( 1200 ID and 
75 km s"' smaller then that of ^ Springob et alj (i2005i) . We attribute 
these differences to our smaller field of view than those of the pre- 
vious works, which have used the Hi 21 cm line integrated within 
large apertures (>3') to derive Vs- The inclination of the disc is in 
good agreement with the one obtained via the expression valid for 
a thin disk: i = acos(b/a) x 44°, where a = 1.3' and b = V are 
the major and minor axes of the galaxy, respectively, as quoted in 
NASA/IPAC Extragalactic Databasjj. The scale length A and the 
bulge mass M are similar to those obtained for the ce ntral region of 
other Seyfert galaxies using similar mo delling (e.g. ISarbosa et al.l 
l2006l : lRiffel&Storchi-Bergmannll20rTI) . 

The stellar velocity dispersion map (top-right panel of Fig.O 
shows a partial ring of low a, values (cr, x 5Q - 60 km s"') at 
r;0'.'6-'.'8 from the nucleus (a230pc) immersed in a background 
of a, X 100 km s"' of the bulge stars. The presence of the partial 
ring is evidenced by an increase in cr, of about 40kms"', from 
locations at 0'.'6 to 1" from the nucleus, while the largest a, un- 
certainties are smaller than 20kms"'. Similar low-cr, structures 
have been observed around the nuclei of other Seyfert galaxies and 
attributed to colder r egions with more recent star formation than 
the underlying bulge JBarbosa et alJl2006l : lRiffel et al.ll2008ll2009l: 



' http://nedwww.ipac.caltech.edu 
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for in te grated spectra liLarkin et alJl l998l: Rodriguez-Ardila et alj 



Figure 11. Rotating disc model for the stellar velocity field (top) and resid- 
ual map between the observed and modelled velocities. The dashed line 
shows the orientation of the line of nodes CFq ~ 112°) and the cross marks 
the position of the nucleus. The color bar shows the range of velocities in 
kms"'. 



iRiffel & Storchi-Bergmannll201 IF) . In particular, a recent study of 
the stellar population in the central region of Mrk 1066 (which 
shows a similar low-cr, ring to that of Mrk 1157) confirms the 
above interpretation, since the stellar population - as derived from 
spectral synthesis - in the low-cr, ring of Mrk 1066 is dom inated by 
intermediate age (10** - 10' yr) stars JRiffel et al1 l2010^. We have 
also performed stellar population synthesis for Mrk 1157, which 
indicates a similar result to that obtained for Mrkl066: the stellar 
population in the low-cr, regions are also dominated by intermedi- 
ate age stars (Riffel et al. 201 1). 

The bulge stellar velocity dispersion can be used to estimate 
the mass of the super-massive black hole in the centre of Mrk 1 1 57: 

log(MB„/Mo) = a +j81og(cr,/o-o), (3) 

where a = 8.13 ± 0. 06, /? = 4.02 ± 0.32 and cto = 200 km s"' 
JTremaine et al.ll2002n . We adopt o", a; 100 ± 10km s"' as rep- 
resentative of the bulge, a value in good agreement with previ- 
ous optical measurements jNelson & Whittle ,1995,) , and obtain 

Mbh = 8.3+^ 2 X lO*^ Mq. 



4.2 Gaseous Excitation 

In order to investigate the gaseous excitation, 
we constructed a spectral-diagnostic diagram 

[Feii]/tl.25yum/Pa/? vs H2/i2. 12/jm/Br7, originally proposed 



|2004 iRodriguez-Ard ila. Riffel & P astorizj |2005[) . This dia- 



gram can be used to distinguish ratios of Seyferts, LINERS 
and Starbursts. In particular, Seyfert nuclei have values in 

the range 0.6<[Fen]/Pa^g2.0 a rid 0.6<H2/Br7<2.0 (e.g. 

iRodriguez-Ardila. Riffel & Pastorizal |2005|). Larger values are 
observed in LINERs and smaller values in Starbursts and Hit 
regions. The diagram is shown in the top panel of Fig.[T2l in 
which the Seyfert-like values are represented by filled circles and 
LINER-like values by crosses. In the bottom panel of Fig.[T2l we 
show how these values are distributed spatially in the galaxy. The 
nucleus has a Seyfert-like value, as well as the emitting gas to 
the east, in regions away from the radio jet. In regions co-spatial 
with the radio jet, the line ratios are LINER-like, and can be 
interpreted as being due to an additional emission of the Ht and 
[Fen] relative to the hydrogen recombination lines due to shocks 
between the radio jet and the gas. This conclusion is supported 
by the [Feii]/Pa/J and H2/Br7 line-ratio maps, as well as by the 
[Pii]/[Feii] map (Fig.|4j, which show the highest values in regions 
co-spatial with the radio structure (see discussion bellow). 

4.2.1 The Hi emission 

The molecular hydrogen excitation can be investigated using its 
flux distribution, line ratios and k inematics and has been s ubject 
of several recent studies (e .g. Storchi-Bergmann et alj 19991: 
Reunanen. Kotilainen & Prietol 2002 : Rodriguez-Ardila et al.l 



2004 [Rodriguez- Ardila. Riffel & PastorizI 



Daviesetai' '2005*; 'Riffel et al.' '2006'. |200"8|, 
Riffel, Storchi-Bergmann & Nagar 2010; Sanchez et al.l 



Ram os Almeida. Perez Garcia & Acosta-Pulidol 12 009: Hic ks et al.l 
i2009l : IStorchi-Bergmann et al.ll2009l l2010l : lOuillard et al.ll2oio[r 
In summary, these studies support that the the H2 emission lines 
can be excited by two mechanisms: (i) fluorescent excitation 
through absorption of soft- UV photons (912-1108 A) i n the 
Lyman and Werner bands JBlack & van Dishoeckl Il987h and 
(ii) collisional excitation due to heating of the gas by shocks, 
due to interaction of a rad io jet with the interstellar medium 
(Hollenbach & McKee 1989) or by X-rays from the central AGN 
( Malonev. Hollenbach & Tielens 1996) . The second mechanism is 
usually referred to as a thermal process since it involves the local 
heating of the emitting gas, while the first one is usually called a 
non- thermal process. 

We can use the H2/i2.2477//i2. 1218 ratio to distinguish be- 
tween thermal and non-thermal processes for the H2 emission of 
Mrk 1157. For fluorescent excitation, typical values for this ra- 
tio are ~ 0.55, while for thermal processes typical values are 
~0.1-0.2(e.g.[M ouri 1994; Reunanen, Kotilainen & Prieto 200j; 
iRodriguez-Ardila et al.l |2004|: IStorchi-Bergmann et all |2009J). As 
observed in Tab.[T] H2A2.2477/i2.1218= 0.11 ± 0.05 for the nu- 
cleus and H2-^2.2477/i2.1218= 0.12 ± 0.03 at a'4 north-west 
from it. Similar values are observed for all locations of the ob- 
served field. This result indicates that the H2 emission is due to 
thermal processes, as observed for other Seyfert galaxies (e.g. 



Rodriguez -Ardila et al.l2004l:lRodriguez-Ardi la. Rjffel & Pastoriza 



2005; Riffel et al 



ioog 



200^ 



201(3 ; IStorchi-Bergmann et alj20 0^ 



Riffel. Storc hi-Bergmann & Nagai 



ermssion can 



The thermal excitation temperature of the H2 
be obtained from the fluxes of its K-band emission lines, through 
the following relation (e.g. Iwilman. Edge & Juhnstone! 20051: 



Storchi-Bergmann et al.l |2009| ; iRiffel. Storchi-Ber gmann & Nagaii 
201(]|) : 
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Figure 12. Top: [Fen]/11.25/jm/Pa(8 vs H2/12. 12//m/Bry line-ratio diagnos- 
tic diagram. The daslied lines delimit regions with ratios typical of Stai'- 
bursts, Seyferts (filled circles) and LINERs (crosses). Bottom: Spatial dis- 
tribution of the emission-line ratios of the top diagram. 



log = constant • 



T, 



(4) 



where F, is the flux of the /''' H2 line, Aj is its wavelength, A, is the 
spontaneous emission coefficient, g, is the statistical weight of the 
upper level of the transition, T, is the energy of the level expressed 
as a temperature and Texc is the excitation temperature. This rela- 
tion is valid for thermal excitation, under the assumption of an or- 
tho:para abundance ratio of 3: 1 . In Fig.[T3]we present the observed 
values for A^upp = j-'- (plus an arbitrary constant) vs £upp = 7", for 
the nuclear spectrum and for the spectrum at position A in Fig.[T] 
(0'.'4 north-west from the nucleus). The fit of the above relation is 
shown in Fig.[T3]as lines and resulted in an excitation temperature 
of T^^c = 2220 ± 50 K for the nucleus and Texc = 2380 ± 45 K for 
the position A. As the observed fluxes are well reproduced by the 
fit of the above equation, it can be concluded that the H2 emitting 



6.0x10^ 




1.2x10^ 



Figure 13. Relation between Wupp = -^ and E^pp = Tj for the H2 emission 
lines for thermal excitation at the nuclear position and at (XM north-west 
from the nucleus. Ortho transitions are shown as filled circles (identified in 
the top fit) and para fi'ansitions as open circles (identified in the bottom fit). 



gas is in thermal equilibrium at Tex^ and that the ortho to para ratio 
is indeed 3, as assumed. 

The result above thus supports thermal excitation for H2. The 
heating can be provided by nuclear X-rays and/or by shocks due to 
the interaction of the radio jet with the interstellar medium. How 
can we distinguish between these two mechanisms? Using the spa- 
tial information and its correlation with the radio: the highest in- 
tensity levels of the H2 emission are elongated following the ori- 
entation of the radio structure (bottom-left panel of Fig.[3]l, what is 
also observed in some of the H2 channel maps of Fig. [TO] Addi- 
tionally, the diagnostic diagram (Fig.ll2t and the H2/Br7 line-ratio 
map (bottom-right panel of Fig.O, show that the highest values for 
YiifBry of up to x 7, are seen in regions which are co-spatial with 
the radio emission, suggesting that the radio jet plays a role in the 
H2 excitation. 

The role of radio jets and X-rays in the H i emission of ac- 
tive galaxies was investigated by lOuillen et al.l (l999 ) using HST 
H2 images of a sample of 10 Seyfert galaxies, as well as their ra- 
dio 6 cm and hard X-ray fluxes. They found no correlation with 
X-rays and a weak correlation with radio 6-cm, suggesting that no 
single mechanism is likely to be responsible for the H2 excitation 
in Seyfert galaxies. This seems to be the case also for Mrk 1157. 
Although we found some spatial correlation between the H2 flux 
and line ratios involving the H2 lines and the radio map, we cannot 
exclude X-ray heating because the H2 emission is less spatially cor- 
related with the radio jet than both the [Fe n] and Pa/3 emitting gas. 
It has also smaller a values (bottom-left panel of Fig.|6j than the 
other emission lines, being similar to the stellar values, suggesting 
that the H2 emitting gas is mainly located at the plane of the galaxy 
and is less affected by the radio jet than the other emission lines. 



4.2.2 The [Fen] emission 

The excitation mechanisms of the [Feu] emission can 
be investigated using the [Fen] /11.2570yum/Paj8 and 
[Feii]/11.2570yum/[Pii]/11.8861yum line-ratio maps shown in 
Fig.|4] The first ratio is controlled by the ratio between the volumes 
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of partially to fully ionised regions, as the [Fen] emission is 
excited in partially ionise d gas regions. In AGN s, such regions can 
be created by X-ray (e.g. ISimpson et alJll996r) and/or shock (e.g. 
iForbes & Wara 19930 heating of the gas. 

For Starburst galaxies, [Fen]/Pa/8< 0.6 and for super- 
novae for which shocks are t he main excitation mecha- 
nism, this ratio is larger than 2 i Rodriguez-Ardila et alj |2004| : 



[Rodriguez- Ardila. Riffel & Pastoriza '2005). As observed in Figs.|4] 
and[T2l for Mrkll57 [Fen]/Pa/3 ranges from 0.6 to 6, with the 
highest values observed at the borders of the radio structure, in- 
dicating that the excitation by radio shocks is important at these 
locations. The stronger correlation between the radio emission and 
the [Fe ii] flux map than with the other emission lines (see Fig.[3l(, 
as well as the higher increase in a at the edge of the radio jet, 
support a higher contribution of the radio jet for the excitation of 
[Fen] than for H2. This result is in good agreement with our pre- 
vious studies of the central region of Seyfert galaxies, which in- 
dicate that the radio shocks may be the main excitation mecha- 
nism of the [Fe 11] near-IR emission lines 1 Storchi-Bergmann et al.l 
I999I: iRilTel et alj|2006l; iRifTel, Storchi-Bergmann & Nagaj|2oiol : 



Storchi-Bergmann et al.ll20 09^. 

The above conclusion is also supported by the 
[Fe ii]/il. 2570 //m/[Pn]/ll. 8861 //m line-ratio map (bottom- 
left panel of Fig.O. These two lines have similar excitation 
temperatures, and their parent ions have similar ionisation poten- 
tials and radiative recombination coefficients. Values larger than 
2 indicate that shocks have passed through the gas destroying the 
dust grains, rele asing the Fe and thus enhancing its abundance and 
thus emission tolivaet alj 2001 : Storchi-Bergmann et al.l l2009l : 
iRiffel. Storchi- Bergmann & Nagan2010n . For supernova remnants, 
where shocks are the domina nt excitation mec hanism [Fen]/[Pn] 
is typically higher than 20 ( Olivaet al.ll200lh . In locations near 
the tips of the radio jet, Mrk 1157 presents [Feii]/[Pii] values of 
up to 10, indicating that shocks are important at these locations in 
agreement with the highest values obtained for the [Fe n]/Pa/S at the 
same locations. In other regions, typical values are [Fen]/[Pn]« 4, 
indicating a lower contribution from the radio jet and higher 
contribution from X-ray heating. 



4.2.3 The coronal-line emission 

Coronal lines are forbidden transitions from highly ionised 
species whose emission extends from the unresolved nucleus 
up to distances between a few tens to a few hundreds 
of parsecs, with the emission-line profiles usually present- 
ing blue wings and being broader than low ionisation lines 
(e.g. 'Rodriguez-Ardila et al\''2006'; 'Storchi-Ber gmann et alj|2009l : 



iMaz zalav. Rodriguez-Ardila & Komossa 2010). 

In the case of Mrk 1157, the coronal-line emission of 
[Caviii] /i2.3220yum and [Si vii] /12.4833;um are spatially resolved 
by NIFS observations, as their flux distributions (bottom-middle 
and bottom-right panels of Fig. |3) present spatial profiles with 
FWHMs;0.'4, which is about 3 times larger than that of the stan- 
dard star. Although the extent of the coronal-line emission region 
of Mrk 11 57 is much smaller than the extent of the low-ionisation 
region, its flux distributions are more extended along the radio jet, 
following the same orientation of the low-ionisation region, sug- 
gesting the same origin. This result supports an origin for the coro- 
nal lines in the inner part of the NLR , as suggested by previous 
authors l lRodriguez-Ardila et al.ll2006l) . This interpretation is also 
supported by the kinematics of the coronal gas, which present sim- 
ilar centroid velocities (Fig.|5j and velocity dispersions (Fig.|6j to 



that of the low-ionisation lines. In addition, the h^ maps for both 
coronal lines show high negative values of up to -0.4, indicating 
the presence of blue wings, in good agreem ent with previous s tud- 
ies of other Seyfert galaxies (e.g. Rodriguez-Ardila et alj2006l) . 



4.3 Gaseous Kinematics 

As observed in Fig.[5]the gaseous velocity field presents blueshifts 
to the east and redshifts to the west, suggestive of rotation, although 
with the orientation of the rotation axis somewhat tilted relative to 
that of the stellar rotation (top-left panel of Fig.|2j. There are also 
differences between the kinematics of the molecular and ionised 
gas, and between that of the ionised H and [Fen] emitting gas. The 
difference relative to the stellar kinematics is largest for [Fen], as 
the orientation of the apparent rotation axis differs from the stellar 
one by more than 45°. 

The gas kinematics and the differences in the velocity fields of 
the different gas species can be better understood by inspection of 
the channel maps of Figs.[8]|9land[T0l Fig.[8] shows that the [Feii] 
emission presents the highest blueshifts to the north-west and not to 
south-east, as expected for the rotation pattern followed by the stars 
(Fig.O. This indicates that these highest velocities are not due to 
rotation. The highest blueshifts, of up to -400 km s"' are observed 
at the location of the 0'.'4 NW radio hot spot. At lower blueshifts, 
there is [Fen] emission also to the south-east, which could thus 
be attributed to rotation in the galaxy plane. At velocities close to 
zero the emission is somewhat more uniformly distributed around 
the nucleus, although it seems to still keep some spatial correlation 
with the radio map. At low redshifts the most extended emission 
drifts to the west, consistent with the rotation contribution, while 
the remaining redshifts to the south-east are observed at the location 
of the SE hot spot, although being spread also to the west of this 
hot spot. At the highest redshifts, there is still emission from the 
structures to the south and from the region of the hot spot at the 
north-west. 

The [Fen] kinematics can be understood as due to the combi- 
nation of the contribution from gas in rotation in the galaxy plane 
and outflows associated with the radio jet, which has the north-west 
part somewhat tilted towards us (in order to produce blueshifts to 
the north-west). At low velocities, the [Fen] kinematics has an im- 
portant contribution from rotation, as evidenced by the emission 
structure moving from south-east at negative velocities to north- 
west of the nucleus at positive velocities, approximately follow- 
ing the orientation of the line of nodes of the stellar velocity field. 
At high velocities the outflows dominate, pushed by the radio jet. 
Nevertheless, the spatial correlation with the radio jet seems to be 
present also in the lower velocity channels, suggesting that there 
is contribution from the radio jet in driving lower velocity kine- 
matics as well. This is also supported by the high velocity disper- 
sion observed over the whole region covered by the [Fe n] emission 
(Fig.lll. 

What is the orientation of the outflow relative to the galaxy 
plane? Inspection of the HST image of Mrk 1157 in combination 
with the observed kinematics, indicates that, if the spiral arms are 
trailing, the near side of the galaxy is the north-north-east (line of 
nodes at PA=112°) and the far side is the south-south-west. The 
near (blueshifted) side of the outflow can thus have three orienta- 
tions relative to the plane: (i) be along the plane, (ii) in front of it or 
(iii) behind it. If it is not along the plane, in order to appear extended 
in the sky, and to show blueshifts to the north-west, it should make 
an angle smaller than 45° with the plane, as this is the inclination 
of the galaxy plane (relative to the plane of the sky), as obtained 
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Figure 14. Re sidual maps betwee n the gaseous and stellar velocity fields. The crosses mark the position of the nucleus, green contours are from the 3.6 cm 
radio image of lNagar et alj 119991) and the color bars show the range of residual velocities in km s"' . 



from the fitting of the stellar velocity field. Inspection of the [Fe n] 
emission distribution at the highest redshifts (350 - 426kms"'in 
the channel maps) reveals also some emission to the north-west at 
the location of the NW hot spot. This suggests that the outflows 
(and the radio jet), although having the north-west side somewhat 
tilted towards us, is launched close to the plane of the sky, so that, 
if the outflow opens in a cone, the blueshifts would come from the 
front wall of the cone and the redshifts from its back wall. If this is 
the case, the cone axis would be behind the plane to the north-west, 
with part of the wall of the cone in front of the plane depending 
on the launching angle and the cone opening angle. In Fig.lISIwe 
present a sketch diagram for the geometry of the outflow relative 
to the plane of the galaxy and plane of the sky. To the south-east, 
the cone axis outflow would be in front of the plane. One difficulty 
with this scenario is the lack of emission to the south-east at the 
highest blueshifts, but this could be due to an asymmetry in the gas 
distribution along the path of the outflow. The hypothesis that the 
radio jet and associated outflow is oriented close to the plane of the 
sky is supported by the Pa/3 channel maps (Fig.|9}, which shows, at 
the highest blueshifts and redshifts, extended emission to both sides 
of the nucleus following the orientation of the radio jet. At lower 
velocities, the main contribution is from gas rotating in the plane, 
from blueshifts to the south-east to redshifts to the north-west. 

The H2 channel maps (Fig.llOl) reach smaller blueshifts but 
are similar to those of Pa/3 for velocities between « -lOOkms"' 
and 100 km s~' , consistent with rotation in the galaxy plane. At the 
highest redshifts there is some emission to the south-south-west, at 
the same location where there is some emission also from [Fe 11] at 
similar redshifts. There is less correlation with the radio structure 
than for Pa/3. 

Distinct kinematics between the H2 and [Fen] kine- 
matics seem to be a common characteristic of active galax- 
ies, having been observed i n man y other Seyfert galax- 
ies JStorchi-Bergmann et alj 19991. Rodriguez-Ardila et al 



2004; [Rodriguez- Ardila. Riffel & Pastorizal 120051; iRiffel et al 



20061 I2OO8I |200^; IStorchi-Bergmann et all I2OIOI : 
Rifi^el & Storchi-BergmannlbOllI) ! 



4.3.1 Comparison between the stellar and gaseous kinematics 

In order to map difi"erences between the stellar and gaseous kine- 
matics, we have constructed the residual maps of Fig.[T4l obtained 
from the difference between the centroid [Feii], Pa/3 and H2 veloc- 
ities and the stellar velocity field. 

The largest and clearest differences are observed in the [Fe 11] 



residuals, which show redshifts of up to 130 km s"' to the south-east 
and blueshifts of up to -80 km s ' to the north-west of the nucleus 
in clear association with the radio jet. 

The PayS residuals are smaller than those of [Fe 11] but there 
are some blueshifts (-30kms"') to the north-east - thus not asso- 
ciated with the radio jet - which appear also in the channel map at 
-43 km s"' (Fig.|9}. The gas excitation is higher at this location (see 
Fig.|4l while the velocity dispersion is lower (Fig.|6l. This could be 
a mild outflow, not shocked by the radio jet and ionised by the 
AGN. There are also some small redshifts (30kms"') to the west- 
south-west, which could be a counterpart of the mild blueshifted 
outflow and some redshifts at the tip of the SE radio hot spot, where 
there is also an enhancement in the gas velocity dispersion, which 
could be part of the the outflow to the SE pushed by the radio jet. 

The H2 residual maps show small redshifts mostly to the 
south-south-west (<40km s"'), a structure also present in the red- 
shifted H2 channel maps (Fig.llOll, which could be part of the mild 
outflow observed in the Pa/3 kinematics, as described above. 

The residual maps support the interpretation discussed above 
for the channel maps. The negative velocities to north-west and 
positive velocities to the south-east, clearly observed in the [Fe n] 
residual map along the radio jet can be interpreted as due to the 
emission of outflowing gas pushed by the radio jet. The [Fe 11] emit- 
ting gas is the best tracer of interaction between ambient gas and 
a radio jet, as the jet destroys dust grains on its way out, releas- 
ing the Iron which produces enhanced [Fe n] emission at the shock 
locations. The presence of shocks explains the high [Fen] velocity 
dispersions. Such enhancements are expected when g as is disturbed 
by a radio jet (e.g. lDopita & Sutherlanj|T99^ . ll996r i. 

The [Fe 11] residual map also supports that the north-west part 
of the radio jet is somewhat tilted towards us, in spite of the fact 
that the channel maps also suggest that the angle between the jet 
and the plane of the sky is small, so that redshifts can also be ob- 
served at the location of the NW hot spot. It can thus be concluded 
that the farthest [Fe n] emission originates in gas which extends to 
high galactic latitudes. This may be the case for some Vap emis- 
sion at the highest velocities which produces some velocity dis- 
persion enhancement seen in Fig.[6] But most of the Pa/J emission 
seems to originate in the galaxy plane, as revealed by the rotation 
pattern which dominates the emission in the Pa/3 channel maps. 
Most of the H2 emission also seems to originate in gas rotating 
in the plane, as supported by the low velocity dispersions (Fig.|6l(. 
The fact that the molecular gas is observed mostly in the galaxy 
plane while the ionised gas is frequently observed to higher galac- 
tic latitudes has also been observed by us for other Seyfert galaxies 
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(e.g.lRiffel et al.ll2006l.l2008ll2009l ; lstorchi-Bergmann et alj|2010l : 
iRifFel & Storchi-Bergmannll201lh . This can be understood as due 
to the fact that, in the plane, the molecular gas is more shielded by 
dust clouds from the ionizing radiation which can destroy the H2 
molecule. 

The velocity fields and a maps for the coronal lines have sim- 
ilar values to those observed for the low ionisation species ([Pii], 
[Fell] and Pa/3), supporting an origin for these lines in the inner 
NLR, in good agreement with the results discussed in the previous 
section. 

In Fig.[T5] we present a sketch diagram for the physical sce- 
nario of the inner 450 pc of Mrk 1157. Two kinematic components 
are observed for the emitting gas. The first is due to gas located in 
the plane of the galaxy presenting a similar rotating pattern to that 
observed for the stars, with a disc inclination of i = 46° relative to 
the plane of the sky and line of nodes oriented along PAs; 112°. The 
second component is due to emission of the outflowing gas within a 
bi-cone, with opening angle of a 60°, oriented along PAr! 153°, be- 
ing oriented approximately perpendicular to the line of sight, mak- 
ing an angle of a 15° with the plane of the sky. 

4.3.2 Mass outflow rate 

In order to quantify the outflows from the nucleus of Mrk 1 157, we 
estimate the ionised-gas mass outflow rate through a circular cross 
section located at 0'.'6 from the nucleus along a bi-cone with the 
same orientation of the radio jet (PA=153°) and opening angle of 
60° (estimated from the left panel of Fig.lI4b. The component of 
the velocity of the outflowing gas along the bi-cone axis (vout) is 
related to the observed velocity (I'obs) by Vom = Vobs/sin^, where Q 
is the angle between the bi-cone axis and the plane of the sky, and 
thus the mass outflow rate can be obtained from 



M,„ 



ImpNeVabJA 



(5) 



sin 6 

where A = nr- a; 3.15 x lO'*' cm^ is the area of the circular cross 
section, nip = 1.67 x lO^^^'g is the mass of a proton, / is the filling 
factor, Nf is the electron density and the factor 2 is included in order 
to account for the outflows to both sides of the nucleus. The filling 
factor can be obtained from the relation 



AnLp^p = 2.046 x 10^'* x N^Vf 



(6) 



where hp^^ is the Pa/3 luminosity of the outflowing gas and V is its 
volume, and the physical quantities are in COS units. We identify 
the outflowing gas as that corresponding to the channels at veloci- 
ties higher than 200kms-'(V >200kms-'and V < -200 kms"'), 
adding the fluxes of the corresponding channel maps plus half of 
the fluxes in the maps centred at -191 kms"' and 181 kms"' (the 
other 50% being attributed to emission from gas in the plane of 
the galaxy), which is 1.3 x 10"'^ ergs"' cm"^, resulting in Lft,^ = 
5.8 X 10'^ erg s ' . The volume of the outflowing gas was estimated 
as a cone with height h = 0.6"= 177.6 pc and opening angle of 60°, 



resulting in V = 5.75 x 10"' cm 



-^ and / : 



: 0.025. 

Assuming N^ = 500 cm"^, / = 0.025 and Voui = 
75 km s"'(from Fig.[l4j, we obtain M^m- 1.55/sin 6 Mq yr"'. 

As discussed above, the radio jet (and thus the bi-cone) is 
probably oriented close to the plane of the sky in order to account 
for the highest blueshifts and redshifts observed in the ionised gas 
emission. If 6 = 15°, Mout ~ 6Moyr"', which is of the order 
of the val ue we have recently obtained f or NGC4151, of M^ut ~ 
2 Mq yr ' JStorchi-Bergmann et alj20IOl). This value is also within 
the range of values reported bv lVeilleux. Cecil & Bland-HawthornI 



(|2005|) for warm ionised gas outflows from luminous active galax- 
ies of Mout -0.1-10 Mayr ' . 

In iBarbosa et al.l ( 12009 ) we have reported much smaller val- 
ues, in the range M^m ~ (1 - 50) x 10"^^ Mgyr"' for a sample of 
six Seyfert galaxies. But in that work we have adopted an electron 
density 5 times lower (A'^ = 100 cm"^) and could not calculate the 
filling factor, adopting a value 25 times lower (lO"""). If we adopt 
the electron density and / values of the present work, we would 
increase the M^ui values by a factor of 150, and the resulting values 
would be in the range M^^t ~ 0.13 - 6.25Mq yr"', more similar to 
the above. 

The same applies to the case of lRifFel et al.l yOOSt), for which 
we used Ng = 100 cm""* and / = 0.01, and the value found of 
Mou, ~ 5 X 10"- Moyr"' should be mul tiplied by 12.5, re sult- 
ing in M„m ~ 0.6 Mq yr"' . In the case of iRifFel et all ( |20IO|) . the 
factor is just 2.5, and the corresponding revised value is Mout ~ 
0.15Moyr"'. 

In the following we compare the mass outflow rate, obtained 
above, with the accretion rate necessary to power the AGN at the 
nucleus of Mrk 1 157, which can be obtained from 



i-bol 



(7) 



where L(,oi is the nuclear bolometric luminosity, 77 is the effi- 
ciency of conversion of the rest mass energy of the accreted ma- 
terial into radiation and c is the light speed. The bolometric lu- 
minosity can be approximated b y L^oi ~ l OOLna, where L hq. is 
the Ha nuclear luminosity (e.g. Ihc1|I999I : IHo et alJ I200II) . The 
Bry nuclear flux measured for a 0'.'25x0'.'25 nuclear aperture is 
Fflry ~ 4.6 X 10"" erg s"'cm"^, as presented in Table[T] As shown 
in Fig.O we obtain a reddening of E{B - V) x 1.7 for the nu- 
cleus of Mrk 1 157 using the Bry/Pa/? li ne ratio. Correcting the Bry 
flux fo r this reddening using the law of lCardelli, Clayton & Mathij 
( 119891) , we obtain Fb,-j, « 8 x 10"'*ergs"'cm"^. Assuming a tem- 
perature T = 10"* K and an electron density rig = 500 cm"^, the 
ratio between Hor and Bry is predicted to be Fua/FBT-y = 103 
dOsterbrock & Ferlanj|20M) . Thus Ln„ ~ 3.7 x 10*" ergs"' for 
a distance to Mrk 1 157 of D = 61. 1 Mpc and the nuclear bolomet- 
ric luminosity is estimated to be Lboi ~ 3.7 x lO"*^ ergs"'. Assum- 
ing T] X 0.1, which is a typical value for a "standard" geometri- 
cally thin, optically thick accretion disc (e.g. lFrank. King & Raind 
l2002h . we obtain a mass inflow rate of m » 6.5 x lO""* Mq yr"'. 

The mass outflow rate (Mout) in the NLR of Mrk 1 157 is about 
four orders of magnitude larger than m, a ratio comparable to those 
observed for other Seyfert galaxies (using the Mout values presented 
above), which indicates that most of the outflowing gas observed in 
the NLR of active galaxies does not originate in the AGN (as the 
nuclear outflow rate cannot be higher than the accretion rate) but 
in the interstellar medium surrounding the galaxy nucleus, which is 
pushed away by the AGN outflow. 

Finally, we can use the above mass outflow rate to estimate 
the kinetic power of the outflow. Followingl Storchi-Bergmann et al.l 
(2010) the kinetic power can be obtained from 



M, 



Ku, + o-~l 



(8) 



where v„„ = v„i,Jsin9 is the velocity of the outflowing gas 
and (T is its velocity dispersion. Using a x 200 km s"' ( from 
Fig.|6j and Vo,„ = v„i,Jsin9 = 75kms"'/>sHil5° = 290kms"' 
(discussed above) we obtain E x 2.3 x 10*' ergs"'. This ki- 
netic power is in good a greement with those derived for the 
Seyfert galaxy NGC4151 dStorchi-Bergmann et aLll2010l) and for 
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Figure 15. Physical scenario of the inner 450 pc of Mrkll57 showing the two kinematic components: rotation in the plane of the galaxy and the outflow, 
represented by a bi-cone. The white-dashed lines show the orientation of the major and minor axes of the galaxy. 



the warm outflow in compact radio sources dHolt et alj201lLl2006l : 
iMorganti et al.ll2005l) . Comparing the kinetic power with the bolo- 
metric luminosity we find that E x 0.15xL(,oi, implying that 15% of 
the mass accretion rate is transformed in kinetic power in the NLR 
outflows. This value is ab out one order of magnitude larger than the 
AGN feedback derived bv lOi Matteo. Springel & Hemquisj J2005h 
in simulations for the co-evolution of black holes and galaxies, in 
order to match the M - a relation and from one to two orders of 
magnitude higher than those derived for six Seyfert galaxies using 
integral field optical spectroscopy (after revi sing the filling facto r 
and gas density values as discussed above) JBarbosa et alj|2009f) . 
Nevertheless, it should be noticed that the E derived here may have 
an uncertainty of up to an order of magnitude due to the uncertainty 
in the geometry of the outflow. 



5 CONCLUSIONS 

We have analysed two-dimensional near-IR J- and A"/ -bands spec- 
tra from the inner » 450 pc radius of the Seyfert 2 galaxy Mrk 1 157 
obtained with the Gemini NIFS instrument at a spatial resolution 
of a535 pc (0'.'12) and velocity resolution of « 35-45 km s"'. We 
have mapped the emission-line flux distributions and ratios, as well 
as the stellar and gaseous kinematics, in the molecular and ionised 
emitting gas. The main conclusions of this work are: 

• The stellar velocity field is well reproduced by a model in 
which the stars follow circular orbits in a Plummer gravitational 
potential, with a major axis at PA=1 12° ± 5°. The near side of the 
galaxy is the north-north-east; 

• The stellar velocity dispersion presents a partial ring (radius 
~250pc) of low values (50-60 km s"') surrounded by higher a, 
values from the bulge stars. This low-cr, ring is interpreted as being 
originated in colder regions with more recent star formation which 



still keep the lower velocity dispersion of the gas from which the 
stars have formed; 

• The stellar velocity dispersion of the bulge (a; 100 kms"') im- 
plies a black hole mass of Mbh = 8.3^2 2 ^ ^^^ ^O' 

• Emission-line flux distributions of molecular hydrogen and 
low-ionisation gas are most extended along PA= 27/153° in agree- 
ment with previous optical [O iii] imaging and following the radio 
jet. This emission extends to at least w450pc (r.'5) from the nu- 
cleus. The coronal lines of [Ca viii] and [Si vii] are resolved extend- 
ing up to a; 150pc (0'.'5) and are slightly more extended also along 
PA= 27/153°; 

• The reddening map obtained via the Pa/J/Bry line ratio 
presents typical values of E(B - V) = 0.5, with the highest val- 
ues of up to E(B - V) = 1.8 seen at the nucleus; 

• The H2 gas is excited by thermal processes and has an exci- 
tation temperature Texc ~ 2300 K. Its emission is mainly due to 
X-ray excitation from the central AGN, but a small contribution 
from shocks produced by the radio jet cannot be discarded; 

• The [Fell] emission shows flux and velocity dispersion en- 
hancements due to shocks produced by the radio jet which destroys 
the dust grains in its way out from the AGN, releasing the Fe; 

• From the comparison between the stellar and gaseous velocity 
fields, as well as from channel maps along the emission line pro- 
files, we conclude that the gas has two kinematic components: (i) 
the first is due to gas located in the galaxy plane, in similar rotation 
to that characterizing the stellar motion; this component dominates 
the H2 emission; (ii) the second component is in outflow, which is 
oriented along the radio jet at PA=153° with the near (blueshifted) 
side to the north-west, but is launched close to the plane of the sky, 
making an angle less than a 30° with the galaxy plane; this compo- 
nent extends to high latitudes and dominates the [Feii] emission; 

• The ionised gas mass outflow rate through a cross section with 
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radius s; 100 pc located at a distance of s; 1 80 pc from the nucleus is 
M„ut~6Moyr-'; 

• Revised mass outflow rates from our previous studies of 8 
other Seyfert galaxies, using the better constrained filling factors 
and gas densities we have determined in the present study, are in 
the range 0.13-6.25 MOyr"', thus similar to the case of Mrk 1 157. 
These rates are much larger than the accretion rates to the AGN, 
implying that most of the NLR outflow mass originates in the in- 
terstellar medium of the galaxy, which is pushed by a less massive 
nuclear outflow. 

• The coronal gas, although being less extended, has similar 
kinematics to that of the low-ionisation gas, suggesting an origin 
in the inner NLR. 

Our final statement, which applies to all Seyfert galaxies we 
have recently observed with NIFS, is that the H2 and the [Fen] 
emitting gas sample different flux distributions and kinematics, in- 
dicating a distinct origin. The molecular gas is in rotation in the 
plane of the galaxy, where is the source of fuel for the super- 
massive black hole. The [Fen] emission, on the other hand, has 
higher velocity dispersions and extends to high latitudes, mapping 
mostly the AGN outflows. 
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